Synthesis of radioactive nuclei and gamma-line radiation 

from novae 



A.D.Kudryashov^, N.N.Chugai^, and A.V.Tutukov^ 
lyiNITI, Moscow, Russia 
^Institute of Astronomy of the RAS, Moscow, Russia 

Abstract 

We carried out kinetic calculations of thermonuclear burning in the 
hydrogen-rich matter to simulate nucleosynthesis yields in nova outbursts. 
These results are used to calculate the light curves of annihilation gamma- 
ray line from N, O and F radioactive isotopes. 

Thermonuclear runaway on accreting white dwarfs (WD) is a conventional 
model of nova outbursts. High temperature hydrogen burning converts most 
of initial CNO nuclei to radioactive isotopes ^^N, i^Q, iSQ, and ^^F (NOF 
isotopes); their presence in novae may be evidenced by 511 keV annihilation 
emission (Clayton & Hoyle 1974; Leising & Clayton 1987). Radioactive ^SNa 
also synthesized in novae may be detected in 1.275 MeV line. Furthermore, 
galactic ^^Al already observed in 1.809 MeV line may be also contributed by 
novae. Here we concentrate only on the synthesis of short-lived NOF isotopes 
and computation of annihilation line flux from novae. 

Following Kudryashov & Tutukov (1995) we calculate nucleosynthesis in 
novae using one-zone model. Temperature and density are assumed constant 
during the burning which is terminated with an exhaustion of the hydrogen 
fraction AX=0.1. The temperature is taken from the range (1 — 3) • 10* K, 
while density is 10"* g cm~'^ in all cases. The adopted composition of CO dwarf 
is X(i2C) = 0.49, X(i60) = 0.49, Xi^^Ne) = 0.01, X{^^Mg) = 0.01, while for 
ONeMg dwarf the composition is X(^^0) = 0.3, X(^°Ne) = 0.5, X^^Mg) = 
0.2. Mixing parameter q (mass fraction of the WD matter in the total envelope 
mass) is varied in the range 0.1-0.9. A kinetic network adopted in this paper is 
an updated version of that from Kudryashov & Tutukov (1995). Apart from H 
and He, it includes nuclei from C to Ca and all nuclear reactions with charged 
particles. 

The obtained amount of NOF isotopes in both types of WD is roughly 
equal to total mass of admixed CO matter of WD in the envelope and thus is 
proportional to the mixing parameter q (Fig. 1). The fraction of ^*F in most 
cases is within the range 10~'^ — 10~^. The mass fraction of ^^Na in ONeMg 
WD envelopes exceeds 1% in a certain region of T — g plane, while in CO dwarf 
envelopes the fraction of ^^Na is usually lower than 10~^. The fraction of ^^Al 
on ONeMg WD for reasonable values of T and q is between 10~^ — 10~^, while 
on CO WD 26 Al content is in the range 10"^ - 10"^. 

We summarize some results of nucleosynthesis computations in Table 1, 
where compositions of major radioactive isotopes are given for average values 
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Table 1: Isotope composition for typical (Tg = 2, q = 0.5, top) and optimistic 
(Tg = 1, q = 0.7, bottom) cases 



WD 


€, 10^' 






ibQ 












erg g-l 


(862 s) 


(102 s) 


(176 s) 


(95 s) 


(158 min) 


(3.75 yr) 


(1.04 • 10'5 yr) 


CO 


4.5 


4.8e-4 


2.1e-l 


2.5e-l 


8.4e-2 


2.4e-3 


2.3e-4 


9.4e-5 


ONcMg 


5.6 


1.6C-4 


l.lc-2 


8.1e-2 


6.0e-2 


1.7C-3 


1.6C-2 


2.0C-3 


CO 


6.1 


l.Oe-1 


2.1e-2 


5.8e-2 


5.4e-3 


7.4e-3 


2.5e-7 


5.5e-6 


ONeMg 


6 


2.3e-2 


4.2e-3 


1.4e-2 


l.Oe-3 


3.9e-2 


2.0e-4 


6.1e-4 



Table 2: Parameters of novae and annihilation line flux (d = 1 kpc) 



Model WD M Vmax /mix -0 F <S> At 
10--'Mq km 10~"t cm-^ s"! cm-^ IQ-* s 



HP 


CO 


2 


2500 


1 


0.31 


19 


53 


2.8 


HHl 


CO 


2 


2500 


1 


0.34 


33 


107 


3.2 


W 


CO 


2 


2500 


1 


0.28 


1 


2.6 


2.6 


HH2 


CO 


2 


2500 


0.34 


0.34 


0.3 


1 


3.7 


HH3 


CO 


2 


2500 


0.99 


0.34 


28 


85 


3.1 


HH4 


CO 


10 


2500 


1 


0.34 


34 


116 


3.4 


HH5 


CO 


2 


3500 


1 


0.38 


70 


214 


3 


HH6 


ONeMg 


2 


2500 


1 


0.25 


17 


56 


3.2 


HH7 


CO 


2 


2500 


1 


0.34 


83 


280 


3.4 


HH8 


ONeMg 


2 


2500 


1 


0.25 


400 


1300 


3.3 



of T and q parameters (typical case) and for parameters, which favour the max- 
imum production of ^^F (optimistic case). Wc suggest that computed isotope 
composition refers to the burning zone, which presumably occupies a fraction ip 
of the envelope mass. This fraction is computed from the energy balance (nu- 
clear energy is the sum of the gravitation binding and kinetic energy) for IMq 
WD. Prior to the ejection, isotopes are presumably mixed in the inner fraction 
/mix of the envelope. 

The flux in the annihilation line for adopted abundances of radioactive iso- 
topes in the typical case (Table 1) is computed assuming that ejecta may expand 
either homologously {v = r/t) or in the form of wind outflow with constant 
velocity and mass-loss rate. The cjccta mass, velocity and /mix are free pa- 
rameters. In the wind case the outflow kinetic luminosity is fixed at Eddington 
limit Lk = (1/2)M?;^ = 10^^ erg s~^. The emergent annihilation luminosity 
is determined by the mass of a transparent outer layer (r < 1). The density 
dependence of a probability of the two-photon positronium annihilation is taken 
into account. 

Given equal outer velocity of cjccta, the emergent luminosity of gamma- 
rays is determined by the density distribution of outer layers. The homologous 
homogeneous sphere (model HHl, see Table 2, Fig. 2a) produces slightly higher 
second maximum (^^F) compared to the homologous envelope with the power 
law density distribution p oc v~'' (model HP). This is caused by the lower 
transparent mass in HP model. The wind model W gives considerably lower 
flux due to the lower density. 
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The effect of variation of mass, velocity, composition (ONeMg WD vs. CO 
WD), and degree of mixing is displayed at Fig. 2b (HH1-HH6 models). The 
minimum mixing degree restricted by the burning zone (/mix = = 0.34) 
severely suppresses both light curve maxima. Remarkably, the marginally in- 
complete mixing, /mix = 0.99, with only 1% of unmixed outer shell results in 
the complete suppression of initial part (t < lO** s) of light curve. Therefore, 
the incomplete mixing, which is quite conceivable in realistic novae makes the 
detection of the first maximum (t < 10^ s) very problematic. 

The Table 2 shows also flux F (distance 1 kpc) in the second maximum 
(t « 2 • lO** s) related to ^*F, integrated flux (fluence $) and the characteristic 
width {At = ^/F). Summing up these numbers, for typical burning case the 
expected fluxes from novae at 1 kpc arc in the range 3 ■ 10"''' — 7 - lO"'' cm^^ c^^. 
The upper limit exceeds by a factor of two the detection limit of BATSE for 0.5 
day-long events (Smith et al. 1996). Given the fact that roughly every three 
year nova may occur at the distance d < 1 kpc the detection of annihilation line 
from novae seems quite plausible. The flux from the optimistic ONeMg case 
(HH8, Table 2), may be as high as 4- 10~^ cm~^ c~^ making such events mostly 
favorable for detection. 
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Figure 1: Radioactive isotope mass fractions in the envelopes of CO WD (top) 
and ONeMg WD (bottom). On the left panel is the temperature dependence 
for gr = 0.5, while shown on the right panel is the mixing parameter dependence 
for Tg = 2. 
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Figure 2: Annihilation line flux from model novae. The left panel (a) shows 
the effect of different kinematics and density distribution for homologous models 
HHl (solid line) and HP (dotted), as well as for the wind outflow model W 
(dashed). The right panel (b) shows the effect of the variation of parameters 
and composition compared to the template model HHl (thick solid line), viz. 
model with no mixing (HH2, short dashes), marginally homogeneous mixing 
(HH3, thin solid line), higher ejecta mass (HH4, long dashes), higher velocity 
(HH5, dash-dotted), ONeMg WD (vs. CO WD) composition (HH6, dotted). 
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